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The concepts of transformation and canonical form have been used in 
analyzing linear systems. In this article we extend these ideas to nonlinear 
systems. A coordinate system and a corresponding canonical form are developed 
for general nonlinear control systems. Their usefulness is demonstrated by 
showing that every feedback linearizable system becomes a system with only 
feedback paths in the canonical form. For control design involving a 
nonlinear system, one approach is to put the system in its canonical form and 
approximate by that part having only feedback paths. 
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I. Introduction 


It is well known that every linear system can be transformed by state 
variable changes and state feedback to a canonical form consisting of multiple 
series of integrators. The number of integrators involved in each series is 
called a Kronecker index, and a linear system is completely characterized by 
its Kronecker indices up to the feedback transformation. We extend this 
result to a large class of nonlinear systems. In a setting of great 
generality, the concept of nonlinear feedback equivalence can be defined as 


follows. A system 

x*f(x,u),xtR n ,ueR m (1) 

is said to be feedback equivalent to another system 

y - g(y,v), y e R n , v e R m (2) 

if there exists a pair of mappings 

T : R n -» R n (3) 

W:R n xR% R ” (4) 


such that (i) T is invertible on a neighborhood near the origin of R n and W is 
invertible with respect to the second variable in R 1 ”, and (ii) for every 
solution (x(t), u(t)) of equation (1), the induced pair of time functions 
(T(x(t)), W(x(t), u(t))) is a solution of equation (2) by substituting T(t) 
for y(t) and W(t) for v(t). The mappings T and W are of course the nonlinear 
generalization of the state variable changes and the state feedback, respec- 
tively, in linear theory. 

According to a result obtained by the authors and G. Meyer [1] (also 
independently by Jakubczyh and Respondeck [2]), a system of the form 

m 

x “ f ( x) + £ g^x) u A (5) 

i-1 


1 


is feedback equivalent to a linear system 


m 

Y - Ay + £ b . v . ( 6 ) 

.,11 
i“l 


if and only if the following conditions are satisfied (with possible 
re-numbering of g . ) : 

1 *r -1 ir.-l 

(a) The set C - {gj.lf.gj] (ad f , g x ), g 2 ,[ f , g 2 (ad 1 f, 

if -1 

g 0 ) , . . . ,g , { f , g J (ad f,g )} spans R , 

4 mm m 

V 2 

(b) The sets C - { g , [ f , g J (ad J f . g. ) ,g» , [f ,g 9 ] , . . . , 

if. -2 3 1 1 if. -2 1 L 

(ad 3 f,g 0 ) g , [f.g ] (ad 3 f,g H are involutive for j ■ 

4mm m 


1,2, • • • , m , 


(c) The span of each Cj is identical with the span of H C, and 

(d) r > r. > . . . > if . 

12 m 

where if, „ , . . . are Kronecker indices for the linear system (6). 

14 m 

Thus, the importance of Kronecker indices has been extended by this 

; 


result to nonlinear systems. Here, however, they play a dual role. In 
addition to controllability span, they also characterize the involutivity of a 
nonlinear system, which has no counterpart in linear systems. 

This result has not only generated a fair amount of interest in the 
analytical problem of feedback linearization (see [1]), but also has been 
applied to the design of nonlinear controllers [3]. A recent successful 
flight test with the UH-1H helicopter by George Mayer at NASA Ames illustrates 
the practicality of our results. In experimenting with such designs, we are 
increasingly concerned with that part of a system's dynamics which is not 
linearizable through feedback transformation. We feel that a measure of the 
"gap" between a nonlinear system and a linear equivalent is needed for 
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robustness analysis of such designs. 


In this article we report our recent progress in this direction. We 
shall present a special coordinate system for the state space of a nonlinear 
system, to be called the s-coordinate system . These coordinates are derived 
on the basis of the intrinsic geometry of a system. In terms of them, a 
canonical form of a nonlinear system is then established. The usual Taylor 
series expansion is replaced by an expansion in terms of the important Lie 
brackets. This form is especially useful for studying feedback lineariza- 
bility. We illustrate this fact by shoving that every feedback linearizable 
system becomes a system with only feedback paths in s-coordinates . We 
emphasize that the s-coordinates can be produced for a general class of 
nonlinear control systems of which the feedback linearizable ones form a small 
(but important) subclass. 

II. s-Coordinate System . 

Let X. ,X 0 ,...,X be an arbitrary sequence of linearly independent vector 
1 i n 

fields in a neighborhood UCR n of the origin. In this section we describe a 
local coordinate system for U which is based on the intrinsic geometry of X,. 
We also assume throughout that the vector fields are analytic. 

We first define our notation and make a simple geometric observation. For 
each vector field X, , let <J> : V -» R n denote the one-parameter group action 

iC K K 

generated by X, , where V is an open set in R x U and includes {0}xU; that 

1C 1C 

is , 

8 

— <j> . (t,p) - X (4>, (t,p)) 

at k k k 

and 
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<*> k (o,p) - p 


for (t.pJeV^. Given any point peU, the mapping <j> k ( • ,p) defines an integral 

curve of X, . More generally, if S is an m-dimensional manifold in U with the 
k m 

property that X. (p) does not belong to the tangent space of S at p for all 
K m 

peS , then the union of all the integral curves $. ( • , p) , peS , defines an 
m K m 

( m+ 1 ) -dimensional manifold S , , and, obviously, X, (q) is a tangent vector of 

m+i k 

s „,i at qcS m +r 

Next we construct a sequence of manifolds associated with X^. We define 


S - {OeR n } (7) 

o 

and Sj_ - {<J> k (t,p) | (t,p)eV k , ptS^^} n U 

where k *> l,2,...,n. From the assumption that X k are linearly independent on 

U and our previous observation, it follow’s that each S, is a k-dimensional 

k 

manifold resulting from the union of the integral curves of Moreover, 

these manifolds form an increasing sequence as 


S C S. C . . . C S C U. 
o 1 n 

We now show that, with the aid of S k> an n-tuple of numbers can be uniquely 

assigned to every point peS^ which turns out to be a local coordinate system 

for S . For a point peS , there is a unique integral curve of X passing 
n n n • 

through p. Suppose this curve intersects S , at a point qeS . . A unique 

n-i n-i 

parameter s n is thus determined by the equation b n (s n ,q)"p. Considering now 


the point qeS^_j , there is an integral curve of X^_j passing through q 


and 


intersection S ^_2 at a point rtS^ ^ . Another unique parameter s n ^ is 
determined by the equation 4>^ ^ ( s n ^ Continuing this process, we 
finally reach the origin of R° with the integral curve of X^ at the n-th step. 
The result is a set of n parameters ( s^ , s^ , . . . , s n ) which is uniquely associ- 


ated with the starting point p and satisfies the equation 

4> ( s n , ( 4> (s 11 1 , ( . . . (<J> , (s 1 ,0) ) . . . ) - p. 

n n— i l 


( 8 ) 
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It is not difficult to see that this mapping F : -> R n , p^ 4 (s^,...,s n ) 

defines an analytic coordinate system for S . This is called the 

n 

s-coordinate system generated by the vector fields 4X. . . 

i- 

In the following, we observe some useful properties of the s-coordinates . 


Lemma 1 . In terms of the s-coordinates, the manifolds defined by (7) are 
simply linear subspaces of R n : 

S k - {s-^.s 2 s n ) e R n |s m -0, k + 1 < m < n} . 

And on each S^, the vector field X^ takes the simple form 

a 

k * 

S, 3s 
k 

or in the form of a column vector 




where X, | denotes the restriction of X, 

N 

This result follows immediately from 
Sj_. For a detailed proof, the reader i< 
noted that the simple expression (9) does 
complete expression for X, is given in thi 


(9) 

on S. . 
k 

the definitions of s-coordinates and 
referred to Spivak [4]. It must be 
not hold for the points not on S^. A 
following Lemma. 


Lemma 2 . In terms of 
expressed as 


s-coordinates and Sj_, the vector field X^ can be 
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* n ( ) * 


X k ' e k * E ^ 


i=l j=k+l i! 


(adV.X, ) 

J K 


>1 


where k * 1,2,. . . , n , and e, - ( 5 S ,) with 6.. being the 

k lk nk ij 

delta . 


Proof. Let X, 


n\ 


\ Kl 


By a Taylor series expansion of X, (s) with respect to s , we have 


1 

fO 

co a i 

K\ 

X k (s) - 

< 

(s 1 S n ,0) + X T 

i-1 i! a(s n ) 

2 

y 

X k 


wj 

! 

IV/ 


/ 1 n 1 <i 

(s , . . . , 9 , 0 j 


From Lemma 1 it is easy to see that the Lie bracket 


ax, ax 

[X ,x, ] - — X X, 

n k ^ n 0 k 

ds ds 


I ax} 


as 


ax: 


k 

as 1 



/ A 


\°l 


( 10 ) 

Kronecker 


(ID 
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3 


We set ( ad°X n> X^) - X^, (ad J X .X^ - [X^.X^J , (ad 2 X .X^ 


(ad 2 X (ad*X , 
n n 


The same computation leads to 


i k 

(ad X ,X.) - . 

n Tc n N i 

3(s ) 


Using the notation of restriction, we rewrite (11) as 

05 (sV 

X. ( s ) * X, I + I (adS: ,X, ) 

k 's . i-1 i! n k S . 

n-1 n-1 


Similarly, X k J can be further expanded with respect to s n * as 


, n-hi 
09 ( S ) 


\| + 1 <adlx n-l'Vl 

's „ i-1 i! 1 


continuing this process until we have 


w ^k+l^i 

\ls -hcls + <adi Vr vl s 

b k+l b k 1 X * b k 


From Lemma 1, X k J - e k> Therefore, combining all the expressions, we obtain 
S k 


n (s J )- 


X,(s) " e v + I X 


i-1 j-k+1 i! 


(adS^.X^ j 
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We have so far only assumed that the are linearly independent on U. 
Additional assumption of involutivenes s will further simplify the expression 
(10). Recall that a set of vector fields is involutive if the Lie bracket of 
any two vector fields is a linear combination of the elements in the set 
with coefficients being scalar functions. 


Lemma 3. If the vector fields X, ,X. ,X are involutive, then X, assumes 

k k+1 n k 

k n T 

the form ■ (0, . . .0,X^, . . . ,X^) in the s-coordinates . 


The proof of this statement is straightforward. We observe that for j ■ 
k+1 , . . . ,n 

(adV.X.) - a°(s)X + c n-1 (s)X . + ...+ a k (s)X. 
j k n n - i k 

from the assumption of involutiveness , where a(s) are scalar functions of s. 

Restricting to we have 


-‘■XL :•**> 


v> 


j-i 


M-i 


M-i 


■ a n e + a 11 ^e , + •••+ ^X . _ 

n n-1 j-1 J-2 


>1 




S i-1 


J- 2 > 


We can further apply Lemma 2 and expand the vector fields or X^^,..., 
to see that these terms must also have vanishing entries above the 


k . 

“ H's. . 

J _1 

k-th place so as to have the form 
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0 \ 


0 

* 


-+ k-th place 


(121 


\*l 

We conclude that each (ad 1 X,,X,)| c has the form (12), and thus Lemma 3 is 

3 R Yl 

proved. o 


If we assume a stronger condition that 

(X.,X k ] - 0 (13) 

for j - k+l,...,n, then X^ - e^. If (13) holds for all l<j,k<n, then a 

classical result follows (see Spivak [4]): X^-e^, k-l,2,...,n. 

In later developments, we also need an expression for a general vector 

field Y on U in terms of the s-coordinate system generated by X ,X_,...,X . 

i / n 

The general expression is given in the next Lemma without proof. 


Lemma 4 . Let Y be an arbitrary vector field on U. In the s-coordinates 
generated by a sequence of vector fields X^ ,X£ , • • • ,X^, Y can be expressed as 


Y(s) 


® n ( s 3 )^ 

Y(0) +11 (ad X X.,Y) 

i-1 j-1 i! 3 



(14) 


IV. Canonical Form 

We now consider the problem of establishing a canonical form for non- 
linear systems 


- f(x) + £ g (x)u , xcR n , 

i-1 


( 15 ) 
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where the vector fields f and g^ are analytic. 

We first construct a sequence of n linearly independent vector fields. 
They in turn generate the s-coordinates on which our canonical form is based. 


We assume that there exists a set of positive integers such that 

l m 

(i) the set C - { gj , (ad*f , gj) , . . . , (ad*" 1 (ad^ *f,g 2 ),..., 

g , ...,(ad m f.g^)} spans R on a neighborhood U near the origin, 

(ii) the span of C ( n C ■ the span of , where C 4 « {g 1 , (ad*f ,g^) , 


j 


j l6 l 

. . . , ( ad^ j 2 f,g.),g,,, .(ad^j 2 f,g_) g , ...,(ad K j 2 f,g )} and 

l z z m m 


j 1 "! * 2 , • • • » , 

( iii) r. >r >. . .>r . 

I Z ID 

Notice that no assumption is made on the involutiveness of the system. 

A crate of m columns is then formed with the vector fields in C as shown 
in Figure 1 . 


g l 

g 2 

. . . 

g m 

(ad*f, gl ) 

(ad*f,g 2 ) s 


(ad*f ,g ) 

* 

• 


- 

(ad‘ r r*f,g 1 ) 





Figure 1 . 


The first column is filled from top down with g^ , ( ad* f , g ^ ) , . . . , 
(ad'l *f,gj), the second column with g^ » ( ad* f , g 2 ) , • • • , ( ad*"2 *f,g 2 ), and so on. 
Now we arrange the vector fields into a sequence in the order from left to 
right and from the bottom row to the top. Hence, we have 
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X 1 - (ad*! l f,g x ) 


X -i 

(a/l'Vgj), 

if 

k 1 >k 2 

X 2 | 

!(ad‘ r 2 _1 f,g 1 ). 

if 

K 1 “ *2 


( ad^3 

if k 

’l" tr 2"‘ r 3 

X 3 " < 

(ad k 'r 2 f, gl ). 

if 

* r l"* r 2>* r 

(ad‘ r 2' 1 f,g 2 ), 

if 



,(adV 3 f, gl ), 

if 

kr l >tr 2 +1 


etc. This procedure can best be illustrated by the following example. 
Suppose we have a system on R^ which forms a crate as shown in Figure 2. 
Then the sequence X^ ,x 2 , . . . .X^ is defined as in Figure 3. 


8 1 

g 2 

g 3 


X 4 

X 5 

X 6 

( ad^f , gj ) 

(ad^ ,g 2 ) 

0 


X 2 

X 3 

0 

(ad 2 f ,gj) 

0 

0 


X 1 

0 

0 


Figure 2. Figure 3. 


Using this set of linearly independent vector fields, we can then 
generate the s-coordinate s . A straightforward application of Lemma 2 and 
Lemma 4 shows the following theorem. 

Theorem 1 ■ In the s -coo rdinates generated by the sequence of vector fields 
previously defined, the system equation (15) assumes the form 


11 



, j \ i 

» n (s J ) 


s - f(0) + X I 


(ad X. ,f) 


i=l j*»l i! 


Vi 


(17) 


n (sV 


+ Z < e k + I I 

k-n-m+1 I i*l j*k+l i! 


(ad 


S'V| 


s j-i 


k-n+m 


where the X^_ are defined in (16). 

We remark that there is an equivalent procedure for introducing the 

i 2 n 

s«(s , s ) coordinates. We simply solve in order the following systems 

of ordinary differential equations with indicated initial conditions. 


dx( s ) 


ds ‘ 


(x 1 , . . . ,x n ) 


. x(0) - 0 


(18) 


, t 1 2. 

dx(s , s ) , 

2 “ X 2 ^ x •••"X ) 

ds 


, x(s 1 ,0) * x( s ^ ) 


. , 1 2 3. 

dxv s ,s ,s ) 


ds' 


X 3 (x J , . . . ,x n ) , x(s 1 ,s 2 ,0) - x(s 1 ,s 2 ) 


1 2 


dx( s ,s ,...,s ) „ , 1 n s , 1 2 / 1 2 n-1, 

■ X (x ,...,x ), x(s s ,...,0) ■ x(s ,s , . . . , s ). 


ds 


If the vector fields X, ,X ot ,..,X are linearly independent, we can solve 

X l n 

X 2 n 

for s ,s ,...,s by the inverse function theorem. 

For our control system (15), the vector fields are defined 

in (16). Since the set C is assumed to span R n , the s coordinates exist for 
our control system. 
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IV. Pure Feedback Systems 


As indicated in a previous discussion, the s coordinates exist under the 
spanning assumptions (i), (ii) and (iii) of Section 3. These assumptions are 
generic, in the sense that "almost all" control systems (15) satisfy them. If 
we also add the involutive conditions mentioned in the introduction, we obtain 
a nonlinear system which is feedback equivalent to a controllable linear 
system. With this extra assumption on involutivity , the s coordinates force 
the nonlinear system (15) into a special form, called a pure feedback system. 
Meyer and Cicolani (5] called this block triangular in a special case and 
showed that it is extremely easy to move from a block triangular system to a 
linear system (the name technique can be used to move from a pure feedback 
system to a linear system). 

In the following definition, we consider system (15) under the spanning 
assumptions (i), (ii) and (iii). We form a crate for the system as described 
in Figure 1. Let n^ be the number of elements in the (»fj-i+l)-th row of the 
crate, and B^-nj+n 2 +. . .+n^, where i,k-l , 2 , . . . . 

Definition 1 . The nonliner system (15) is a pure feedback system if it is 
of the form 


•1 

x 


- 1 , 1 

f (s 



) 


•2 

x 


1 2 , 

f (x X ) 


• B 1 B l, 1 B 2, 

c « f (x,...,x ) 


• V 1 A* 1 , i \ 

X - f (x , . . . ,x ) 


13 




-\ 


,B,+2 B +2 B, 

X - f (x .... ,x ) 


.^2 1 ^3 

x - f (x x ) 


B 0 +l B 0 +l , 

• k.- 2 r.,-2 ,1 n, 

x 1 - f 1 ( x , . . . ,x ) 


• n— m -.n-m, 1 n, 

x ■ f (x , ... ,x ) 


•n-m+1 

x 


..n-m, 1 
f (x . 


. ,x ) 


m -i 

, n-m+1., 1 

+ I -84 < x 
i-1 


n, 

»x )u i 


m 


»* t * > T* J 

X - f (x , . . . ,x ) + l g . (x 

i-1 


n \ 

.X ^ 


We remark that B -n-m, B -n,+...+n -n, and n <n <n < . . . <n 

rj-1 k - ! 1 r l 12 3 k-j 

definition is motivated by the following diagram. 




W”/ 



l x 


Figure 4. 


This 


\ 

7 
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We note that there Is no feedforward signal flow other than the main 
path. Consider again a system with an associated crate as shown in Figure 2. 
It is a pure feedback system if it is of the form 


f (x^x^.x^) 


f (x^ , Sg ' • • • .x ) 


f 3 (Xj ,x^ , . . . ,x^) 


A, 1 6. i 4, 1 

f (x , . . . ,x ) + l g , Cx , 

i-1 


6 , 

,X )u_^ 


- 5 , 1 6 \ J. V 5 / 1 e \ 

f (x X ) + l g (x , . . . ,X )u . 

• 1 1 1 
1-1 


, 6 , 1 6 , ~ 6 , 1 
f (x X ) + X g,( x ’ 

i-1 


6 * 

.X )u. 


To the spanning assumptions (i), (ii) and (iii), we add 

r .-2 

(iv) The sets C - {g ,[f,g ] (ad 3 f , g, ) , g 9 , [ f , g 9 ] , . . . , 

K -2 3 1 r .-2 

(ad f,g„),...,g ,(f,g ],...,( ad 3 f,g H are involutive for 
l m m m 

u-1 , 2 , . . . ,m. 

Under this additional assumption (iv), system (15) is feedback lineariza- 
ble [1]. We show that it takes a special form in the s-coordinates . 


Theorem 2 . If the sets of vector fields C ^ , j-l , 2 , . . . ,m are involutive, then 
the nonlinear system (15) is a pure feedback system in the s-coordinates. 


Proof. For our nonlinear system, the vector fields X, ,X„,...,X are defined 
i l n 

from those elements in the set C by Equation (16). We take as before 
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rij “ number of elements in C with ad superscript ir^-l, 
n 9 ■ number of elements in C with ad superscript r^-2. 


n - m ** number of elements in C with ad superscript 0. 

K l 

Also, we have - n^+n^t . . . +n^ , l<k<r^ with B^ “ n. The canonical form 
for the system in the s-coordinates generated by is 


Jv 1 


n ® ( s J ) 

f(0) + l l ( ad i X . , f ) 


j-1 i-1 i! 


n 

+ I 

k-n-m+1 


(17) 


1-1 

n 00 (sV 


♦ I I 

j-k+1 i-1 i! 


(adVj ,X k ) 


k-m+m 


’i-1 


It has been shown 
that each of the sets 


in'|lj that the involutivity of C_.,j-1,2, 
of vector fields {X +1 X ],k-l,2 


,m implies 
. , r j-1 is 


involutive. 

In particular, when k-K^-l 
involutive. From Lemma 3, this 
the form 


, the set {X + 1 , . . . ,X } - {g,,...,g ] is 
p, m 1 m 

k 

implies that the control terms in (17) is of 


1°) 


1 0 ^ 


l°\ 

0 


0 


0 

* 
i : 

u i + 

* 

U 2 + • • • + 

* 

\* 


\*i 


{* 


( 20 ) 


where * means possibly nonvanishing functions of s ^ , S£ , • • • . s^, and they begin 
at (B ,+l)-th entry. 

V 1 

The second summation is rewritten as - 
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Bj (s 3 ) 1 

I I 

3-1 i-1 i’- 


(ad X X. ,f) 


B 2 ® (s' 1 ) 1 

+ I I 

j-B-^+1 i-1 i! 


(ad 3 *. ,f) 


+ l I 

j-B i+l i-1 i! 
J r — 1 

i 


(adV ,f) 


Observe that as the index j runs from B , +1 to B -n in the last term 

V 1 J r i 

of the summation (21), the vector fields (ad X^,f) are either simply 

X. ,6 _+l<k<B or linear dependents on X.,B „ + l<)i.<B • Further, because 

k JT ^ Z ^ ^ Z ^ J. 

the set -{X X -X } is involutive, it is easy to see that 

p n + 1 p n 

V 2 


(ad 1 X^ , f ) - 



+- (B „+l)-th place, 

V" 


When restricted to the asterisks in (22) are functions of 

s^,...,s n , As we sum up the series, we have 

\ • i i 

I I (ad x x f) 

j-B ,+l i-1 i! J ‘S. . 

J tfj-1 3-1 


+- (B _« + l)-th place 
*1 
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where * are functions of , s ' , . . . , s 11 . 


Continuing the analysis with every term in summation (21), it is seen 


that 


(ad 1 X^,f) - 


0\ 


I*/ 


<■ (B , +l)-th place 

K* ^ Jv 


(23) 


when the term with j running from B , +1 to B is considered 

j ™ j ” K ' 1 

S. , is the manifold 

J-l 

sterisks in (23), when 
^ ^ . Summing up the 


where k“l , 2 , . 

. . , k j , and 

B 0 -0. In 

thi s 

case , 

characterized 

by s^“s^ + ^- 

n n 

. . . ■ s *0 . 

Thus , 

the a 

restricted to 

S . . , are 

functions 

, 1 
of s 

i • • • » s 


subseries, we have 

B , +1 ® 


-k 

X 1 


(.b 1 , 

I ( ad a , f ) 

i J 


i-1 


1 ! 


S . 
J 


I o\ 


<- (B ,+l)-th place 


*/ 


B , +1 

1 2 Kj-k ^ 

where * are functions of s ,s ,...,s , but is independent of s , 


"• 

This concludes our claim that the system is a pure feedback system, o 


We remark that for this paper we have worked in a neighborhood of the 
origin in R n , but any point in R n where the system satisfies the spanning 
conditions could have been considered. 
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V. 


Conclusions and Related Research 


For a real analytic system of the form 
m 

x - f(x) + Y g.(x)u. , xtR n 

i=l 1 1 

we have introduced a coordinate system which reflects the intrinsic differen- 
tial geometry of the system. The basic vectors of the s-coordinate systems 

are generated by the vector fields in the set C, and the indices r. ,ir 

i i m 

correspond to Kronecker indices for a controllable linear system. Moreover, 
we have introduced a "canonical expansion" of the nonlinear system in the 
s-coordinates, and the usual Taylor series terms are replaced by a Lie bracket 
expansion. Under involutive assumptions on subsets of the vector fields in C, 
we show the nonlinear system is a pure feedback system in the s-coordinates. 

In general, a "measure of how close" a nonlinear system is to a pure 
feedback system should be indicated in its canonical expansion in the 
s-coordinates. This is especially interesting since pure feedback systems are 
feedback equivalent to controllable linear systems. 

It is known that in the s-coordinates the usual Taylor series linear 
approximation to a nonlinear system and a linear approximation generated by 
certain Lie brackets agree [6]. Adding outputs to a nonlinear system allows 
us to discuss Volterra series expansions. Since Volterra series kernels 
involve Lie derivatives and Lie brackets [7], the s-coordinates may be natural 
for such series. Moreover, the s-coordinates may prove advantageous for 
nonlinear input-output systems realization. 
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